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Abstract We attempt to calculate the gravitational time delay in a time-dependent gravi- 
tational field, especially in McVittie spacetime, which can be considered as the spacetime 
around a gravitating body such as the Sun, embedded in the FLRW (Friedmann-Lemaitre- 
Robertson- Walker) cosmological background metric. To this end, we adopt the time transfer 
function method proposed by Le Poncin-Lafitte et al. (Class. Quant. Grav. 21:4463, 2004) 
and Teyssandier and Le Poncin-Lafitte (Class. Quant. Grav. 25:145020, 2008), which is 
originally related to Synge's world function Q{xa,xb) and enables to circumvent the inte- 
gration of the null geodesic equation. We re-examine the global cosmological effect on light 
propagation in the solar system. The round-trip time of a light ray/signal is given by the 
functions of not only the spacial coordinates but also the emission time or reception time 
of light ray/signal, which characterize the time-dependency of solutions. We also apply the 
obtained results to the secular increase in the astronomical unit, reported by Krasinsky and 
Brumberg (Celest. Mech. Dyn. Astron. 90:267, 2004), and we show that the leading order 
terms of the time-dependent component due to cosmological expansion is 9 orders of mag- 
nitude smaller than the observed value of dAU/dt, i.e., 15 ±4 [m/century]. Therefore, it is 
not possible to explain the secular increase in the astronomical unit in terms of cosmological 
expansion. 
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1 Introduction 

Technological progress has significantly enhanced the accuracy of solar system observa- 
tions. Presently, to archive an accuracy of 9 to 11 digits or higher, the 1st order ff{c^^) or 
^(G) (c is the speed of light in vacuum and G is the gravitational constant) post-Newtonian/post- 
Minkowskian effects of general relativity must be incorporated in a number of observational 
models. Moreover, several space missions such as GAIA [1|, Laser Astrometric Test Of 
Relativity (LATOR) (2), and Astrodynamical Space Test of Relativity using Optical De- 
vices (ASTROD) (21 are planed and especially LATOR and ASTROD can potentially detect 
post-linear or 2nd order i^(c^^) or GifP') post-Newtonian/post-Minkowskian effects. 

Highly accurate and rigorous light propagation models are required to analyze high 
precision data such as the radar time delay, deflection of light rays, and frequency shift 
of radio signals. Although, these effects have been extensively investigated (see i.e. (4]|5l|6] 
|7]l8l|9l and references therein, most of the previous studies were based on integrating the 
null geodesic equation of light rays. This approach is effective in the case of linear or 1st 
order approximation. However, it is generally difficult to compute the contribution of 2nd 
and higher order effects, even in the case of static gravitational fields; the situation becomes 
further complicated in the case of time-dependent spacetime. 

A novel method known as the time transfer function method has been recently developed 
lilOllllI to overcome such circumstances. This approach is originally based on Synge's world 
function 12 (;cyi,jCB) I12II . 



where g^v is a metric tensor of spacetime, = (x^ = c?^,x^ = x^) and xb = (jc^ = crg,;c^ = 
xg) are the coordinates of end points A and S, respectively, on the geodesic world-line, and A 
is the affine parameter. Then, the world function 12 (xf^^x^) is defined as the half length of the 
world line between A and B. On the basis of ([T]l, the time transfer functions introduced in IIIOI 
[TTl define the travel time of a light rays/signal between two points, i.e., 1^ — 1/^. \W\ and 1111 
succeeded in constructing the general A'-th order post-Minkowskian approximation scheme 
with respect to the gravitational constant G via an iterative procedure, and they showed 
that this method can reduce the number of calculations by taking the static Schwarzschild 
spacetime for instance. 

On the other hand, in the field of experimental relativity, astrometric observations in 
the solar system have played a crucial role in the verification of gravity. It has been show 
that the main PPN (parametrized post-Newtonian) parameters jS and 7 are equal to unity for 




(1) 
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general relativity, (jS = 7 = 1); as increasing accuracy, jS and y are strictly constrained to 
unity (see Fig. 1 of (13\ ). The enhanced accuracy of observations accounts for the existence 
of unexplained phenomena by means of the framework of general relativity; such phenom- 
ena include the anomalous acceleration of the Pioneer spacecraft 1141 . Earth flyby anomaly 
1151 , the anomalous perihelion precession of Saturn |jT6|, and the secular increase in the 
astronomical unit 1 17 ). Presently, the origins of these anomalies are far from clear; however, 
they may be attributable to some fundamental properties of gravitation (see I18II19I and the 
references therein). 

The secular increase in the astronomical unit is particularly interesting. Here, the astro- 
nomical unit implies that the conversion constant between two length units and the current 
best-fit value is 1201 

= AU = 1 .495978706960 ±0.1. (2) 

1 [mj 

The determination of astronomical unit is directly related to the light propagation formula, 
i.e., the equation of light time. Without loss of generality, the secular increase in the astro- 
nomical unit was found by the relation 



'theo — 



dAU , 
AU+-j-(f-fo) 



[s], (3) 



where ftheo is the computed round-trip time (theoretical value), d^eo is the interplanetary 
distance obtained from planetary ephemerides (note that the length unit of t/theo is [AU]), 
AU and dAV/dt are the astronomical unit and its time variation, respectively, and to is 
the initial epoch. f,heo is compared with the observed round-trip time fobs- Krasinsky and 
Brumberg estimated dAl]/dt as 15 ±4 [m/cylQ. 

The influence of cosmological expansion has been investigated as a potential cause of 
secular variation in AU Ill7ll21ll22l . In particular, 11171 and f22\ examined its contribution 
to light propagation. The approach of 1171 is somewhat qualitative, and the discussion of 
11221 is essentially reduced to a static model that is equivalent to Schwarzschild-de Sitter 
spacetime. These circumstances are due to the fact that the time-dependent null geodesic 
equation must be solved in order to study the cosmological effect on light propagation; 
however it is generally difficult to solve this equation. Nonetheless, the time transfer function 
method proposed in IIOIIIII is applicable not only to the static case but also to the time- 
dependent case. Therefore, it is useful to re-examine the effect of cosmological expansion 
on the basis of the time transfer function and to compare the theoretically obtained results 
with the observed dAV/dt. 



In tiiis paper, cy denotes century according to Krasinsky and Brumberg (2004). 
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In this paper, according to the time transfer function method, we consider light propaga- 
tion in a time-dependent gravitational field, and we focus our discussion on the gravitational 
time delay. We adopt the McVittie model I231I24I as the time-dependent spacetime; it can 
be considered as the spacetime around a gravitating body such as the Sun, embedded in 
the FLRW (Friedmann-Lemaitre-Robertson- Walker) cosmological background metric. The 
time-dependency in the solution is characterized by either the emission time or the reception 
time of a light ray/signal. Then, the obtained results is applied to the secular (time) variation 
in the astronomical unit, reported in 11171 . 

The remainder of this paper is organized as follows. In Section[2l we briefly summarize 
the time transfer function, and in Section |3] we investigate the gravitational time delay in 
McVittie spacetime. In Section|4l we discuss the secular increase in the astronomical unit. 
Finally, in Section[5] we conclude the paper with a summary of our study. 

2 Time Transfer Function 

Before discussing the gravitational time delay in McVittie spacetime, we shall briefly sum- 
marize the time transfer function method (see lllOlllll for the details). 

It is generally difficult to acquire the world function Q{xatXb) concretely; nevertheless, 
in the case of Minkowskian flat spacetime, the world function is easily obtained by using 
the parameter equation x{X) = {xg —xa)^ +xa and by setting = and = 1 | [10lll2ll : 

n^"HxA,XB) = \iM4-4)i4-4)^ (4) 

where the jr^ (jU = 0, 1 , 2, 3) are Minkowskian coordinates with respect to the Minkowski 
metric r)^v = diag(— 1; 1; 1; 1). For the null geodesic, the world function Q (xa,xb) satisfies 
the condition 

n{xA,XB)=0 (5) 

because ds^ = 0. Hence, from ^ and ([S}, the travel time between A and B, i.e., — Ia, in 
Minkowski spacetime becomes 

c^itB - tA f = 5,v(4 -4) = I^Ib^ (6) 

where 5,^ is Kronecker's delta. The time transfer function approach starts from (|6}, and the 
weak-field approximation is recursively developed with respect to the gravitational constant 
G. 

If the metric has the form 

gjiv = 'r\iiv+hiiv, |/!mv|<1> (7) 
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where /i^y is a perturbation to r]^v, the time transfer functions that give the travel time of 
the light ray/signal are formally expressed as follows: 

tB-tA = ^e(f/i,XA,XB) = -[RAB+\{tA,'^A,'X.B)] (8) 
C 

= .9'r(xA,tB,'!i-B) = - [RaB + M^aJb ,'a-B)] , (9) 

c 

where ^(/^,Xyi,XB) is the emission time transfer function in spacial coordinates xa,xb and 
tA is the emission time of the signal, ^(x^jf^jX^) is the reception time transfer function in 
spacial coordinates xa,xb and Ib is the reception time of the signal, Rab = \xb — xa\, and A^. 
and Ar are the emission time delay function and reception time delay function, respectively. 
Ae and A^ characterize the gravitational time delay. In ([8) and (|9), Rab comes from (|6}. 
Henceforth, A denotes the emission and B denotes the reception. 

Ae and A^ can be iteratively calculated up to any order of approximation IIOIIIII : never- 
theless, in this paper, we only need the 1st order formulae, i.e., 

= ^ [sw - ^N^Bsli) + A'AXb^;{)] dn (10) 
= ^ /o' [8fi)-2NU^i^+N!,BNiB8[{]dX, (11) 

where indicates the 1 st order perturbation with respect to r; and A'^^ = {x'^ — x\ ) /Rab ■ 
In jlOb and jllb . integration is carried out along the straight line 

t{^)=tA + ^iTAB, x{pl)=xa+Pl{xb-xa) for Zii'', (12) 
t{X)=tB-XTAB, x{X)=xb-)i{xb-xa) for Zii'', (13) 

where Tab is the time lapse between A and B along the straight line. Then, we can put 
Tab = Rab/c- 



3 Gravitational Time Delay in McVittie Spacetime 



3.1 McVittie Spacetime 



The McVittie metric is expressed in standard comoving form 02311241 as 



1 



2c^ra(/) 



1 + : 



C^dt^ + 



1 + 



GM 



2c^ra(t] 



a^{t){dr^ + r^dQ^), 



(14) 



where dQ^ = dG^ + sin^ 6d^^, M is the mass of the central gravitating body, and a{t) is a 
scale factor. l ll4b reduces to the Schwarzschild solution when a(t) = constant, and it reduces 
to the FLRW cosmological model for the curvature parameter k = when M = 0. 
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The various observational models of the solar system are currently formulated in some 
kind of proper coordinate system such as the barycentric celestial reference system (BCRS) 
based on the post- Newtonian framework 1251 , instead of the cosmological comoving frame. 
Hence, to compare the effects formulated using the proper coordinates with the cosmological 
ones within the same framework, we adopt the radial transformation 

f rlUI ^ 

R = a{t)r 



GM 

1 + - 



2c2ra(r) 



(15) 



then l ll4t is rewritten as 




/ , 2GM\ 2,2 dR HR \ ,^2 
ds = - \\ ^ &-dt^ + , cdt +R^dQ^ 



1 + ^ + ^(M^) ) cdtdR 



2GM 



1 + + ff{M^) I dR^ +R^dQ^, (16) 
c^R 

where H = H{t) = d{t) /a{t) is the Hubble parameter. 

As mentioned in previous section, the light path used in the computation is rectilinear 
so that the rectangular coordinate system can be used instead of the spherical coordinate 
system. By coordinate transformation, 

X = 7?sin9cos(^, >' = /? sin9sin(^, z = /fcos0, (17) 

and I ll6b becomes (see |31 1) 

, / 2GM H^R^\ J J 2Hx' ( 2GM ,A 

2GM ■ ■ T \ . . 

5,/ + -^^x' + eiM^) dx'dx', (18) 
c-^R^ I 



where R = \J x^ +z^. To simplify the computation, the straight line used in integration 
is parallel to the ;i;-axis (see Fig.[T}. Hence, ilSi reduces to 

, / 2GM H^R^\ , , 2Hx f 2GM ^„ ,A 
ds- = -(l--^--^)c\U^-—(x + -^ + m'))cdtdx 

+ fl + ^,2 + ^(M2)Vx^ (19) 



where y = b = constant {b is the impact factor), z = 0, /? = \/x^ + b'^, and in this case, we 
may put N^^ = {xb — xa) / Rab , N^ab ~ ^ab ~ ^- "^^^^ approach is similar to the one described 
in 1321, (see Section 40.4 and Fig. 40.3 of (32)). In we see the Hubble parameter H{t), 
which generally is an arbitrary function of time t. Hence, we suppose that H{t) changes 
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Actual light path 



Light path used 
in calculation 




Fig. 1 Light ray/signal patli. We assume tliat tlie patli used in calculation is rectilinear and parallel to the 



X-axis (bold solid line); b is the impact parameter and R = Vx^ + b^. The actual path is denoted by the bold 
dashed line. 



adiabatically because the time lapse of a light ray/signal observed in the solar system is 
much shorter than the age of the universe, Tu w 10^" [yr]: 



Hit) =Ho + 



dH 

IF 



t + > 



cfH 



dH 



dt 



Ho 



10 



-24 



[l/(syr)], 



(20) 



where //o « 10"'^ [1/s]. 



3.2 Gravitational Time Delay 



The time delay functions A^. and A,- can be resolved into the following components: 

A,{xA,tA,XB) = A{xA,XB)+A{xA,XB)+A,{xA,tA,XB) + G{M^,{H\Q)^,H\Q) , (21) 
ArixA,XB,tB) = A{xa,Xb) + A{xa,Xb) + Ar{xA,XB,tB) + G {M^ ,{H\of -.Hlo) , (22) 

where A{xa,xb) corresponds to the Shapiro time delay due to central mass M, A{xa,xb) is 
due to the static component of cosmological expansion (Hq in il20b). and Ae{xATtA,XB),Ar{xA,XB,tB) 
are due to the time-dependent (secular) component of cosmological expansion (dH/dt\Q in 
ilOi). By the straightforward calculations from (llOt . (1111 . i ll9b . and ( 120b . we obtain follow- 
ing results: 



A[xa,xb) = —r 



2 In 



XB+Jxl+b^ 



XA+ \/x\+b^ 



xb 



xa 



/xl + b^ Jxl + b^ 



(23) 



c 

4GMHo 



A{xa,xb) = -e-^{xl-xl) + -^[xl-xl + 3b^{xB-XA)] 



xl + b^-Jxl+b^]+G {m\ {H\o)\H\o) (24) 
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Ae{xA,tA,XB) 



{xB-XA}tA + :^{2xb - xaxb - xa)Tab 



dH 
dt I 

+ ^{[xl-x\+'ib'^{xB-XA)]tA 

+ 7 [3-^1 -xa^^b- XaXb -xl + 6b^ {xb -xa)] Tab } 



[{xb - 2xa)Tab + 2{xb - XA)tA\ \/xl + b^ 



2GM 

c^Rab 



+ [xaTab -2{xb- XA)tA] \Jxl+b^ 
-b^TABln 



Xb + \/xl + b^ 



XA + \/xl+b^ 



+ ff{M\{H\of,H\o) 



X , , dH 

Ar(XA,XB,tB) = — 



1 



{xb -xD'b- -^{xl+ xaxb - 2x1 ) '^AB 



+:^{[xl-x\+3b'^{xB-XA)\ tB 



- - [x| + xax^b +x\xb- 3x| + 6b^ {xb -xa)\Tab] 



2GM 



+e^- <^ \xBTAB-2{xB-XA)tB] Jxl+b"^ 



■ [{2xb-xa)Tab - 2{xB-XA)tB] Jxl+b^ 



-b^TABln 



XB + \/x^g + b2 



XA + \/x\ + b^ 



+ff{M\{H\o)\H\o): 



(25) 



(26) 



where e = A';^^ = 1 for xb — xa > and e = — 1 for jc^ — -^a < 0- £ is derived from the term 
-2NXBg°[-^ in Onil and O- 



3.3 Equation of Light Time 

Actually, xa and xb are positions at r = and t = ts, respectively: 

XA=XA{tA), XB=XB{tB)- (27) 

Then, ll21t and ll22t can be symbolically expressed as 

Ae{xA{tA),tA,XB{tB)) = A{xA(tA) ,XB{tB)) + A{xA{tA) ,XB{tB)) 

+A,(xA{tA),tA,XB{tB)), (28) 

^r{xA{tA),XB{tB):tB) = 4 (x/i (f^) , ) + 4 (x^ (Ja) , ) 

+ArixA{tA),XB{tB),tB), (29) 
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(-'(2) 



*-'(q) 



t' 



(N) 



Reflector 



Emitter/ 




Time >■ 



Receiver ^(i) ^"'^'f'^) ^''^^'t'"*' t"'''' t"^' 



'(N) 



Fig. 2 Time-dependency of solutions. The emission time, reflection time, and reception time of the (</)-th 
observation are denoted by and/"'''), respectively. 



Hence, the equations 



tB-lA = - [RAB{xA{tA),XB{tB)) + \{xA{tA),tA,XB{tB))], 
C 

tB-tA = -[RAB{xA{tA),XB{tB))+ K{xA{tA),XB{tB),tB)\, 



(30) 
(31) 



can be considered as the equation of light time. Here that xa and must be known before- 
hand as a function of time t via the numerical integration of the equation of motion or via 
lunar-planetary ephemerides. 



3.4 Time-dependency of Solution 

The solutions from (121b to ( I26t have the following form: 

,1^(9)^ J?) f(i) _ 1 \j)(i)(Si) M)\<Ai)(M) M)\ 

\Xj^ ,Xg j — ^ Y^ABy^A '-^B l^'^l V^A '"^B ) 

+E^^\x''^\x^^)t^^], (32) 

\Xj^ ,Xg ,lg ) — ^ ["aB^.-^A '^B \Xj^ ,Xg ) 

+Ri'\x'f\xi'^)ti'^], (33) 

where the superscript (q) implies the cj-th observation, E\''\xf\x'g'^),R^^\x^^\x^^'') consti- 
tute the static component of gravitational time delay, and E2''\x^^\x^^'')tA,R^2'' {x''^\x''g'^)tB 
constitute the component depending on either or fg. If we carry out a series of observa- 
tions, q = 1,2,3, - ■■ ,N (see Fig[2ll, /^\/g ^ on the right-hand side of ( 132b and ( 133b can be 
regarded as the time-dependent (secular) component of the solution. 
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4 Application to Secular Increase in Astronomical Unit 

The Secular increase in the astronomical unit (AU), reported by Krasinsky and Brumberg 
(2004) [17], is an unexplained physical phenomenon observed in the solar system. This 
anomaly was discovered while analyzing planetary radar and spacecraft (mainly Martian 
landers/orbiters) ranging data and improving the various astronomical constants including 
AU. [nl estimated 



as the most appropriate value. Subsequently dAV/dt ~ 20 [m/cy] was separately evaluated 
by Pitjeva at the Institute of Applied Astronomy (lAA), Russia, and by Standish at the Jet 
Propulsion Laboratory (JPL), USaB 

Previously, some attempts have been made to explain the secular trend in terms of cos- 
mological expansion I17ll21ll34l . In particular, [1 17 11221 considered its contribution to light 
propagation. However, as mentioned in Section [T] it is generally difficult to compute the 
time-dependent geodesic of null rays; hence, the approach of the former is somewhat quali- 
tative, whereas that of the latter is essentially restricted to discussion in static spacetime. 

Now, on the basis of the results obtained in previous section, let us re-examine whether 
the cosmological effect relates with the observed dAU/dt. It is appropriate to regard the 
coefficients of and Ib in ( 132b and ( 133b as secular terms owing to cosmological expan- 
sion. Then, if we assume that dH/dt\() w 10^^^ yr)] from ( 120b , the leading order of 
magnitude of coefficients is. 



because xa and are of the order of a few [AU] or 10^' [m] in the solar system. Unfor- 
tunately, this is approximately 9 orders of magnitude smaller than the evaluated dASJ jdt, 
i.e., 15 [m/cy]. Therefore, the time-dependent effect due to cosmological expansion does 
not induce the secular increase in AU. 

5 Summary 

We considered light propagation in the time-dependent McVittie spacetime, which can be 
considered as the spacetime around a gravitating body such as the Sun, embedded in the 

^ 1171 used only the radiometric obsereational data, whereas Pitjeva and Standish carried out analysis using 
various kinds of data such as optical and VLBI; each of them worked with different software. However, as 
it stands, their results have not been published officially because of the difficulty in interpreting the obtained 
value, i.e., 20 [m/cy] f331. 




(34) 



£2(-«a,^b)w^2(-«a,xb) wlO ^" [m/yr] = 10 [m/cy] 



(35) 
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FLRW cosmological background metiic. To discuss the time-dependent null ray, we adopted 
the recently developed time transfer function method lllOlll II which is originally related to 
Synge's world function Q{xa,xb) and precludes the integration of the null geodesic equa- 
tion. The first step involved the application of this method to McVittie spacetime and re- 
examination of the cosmological effect on the round- trip time of a light ray/signal in the 
solar system. The time delay functions Ag and Ar, which characterize the gravitational time 
delay, were given by the functions of not only the spacial coordinates xa and xb but also the 
emission time or reception time (see ll21l to i26i): the presence of the terms and Ib in 
zip and A,- express the time-dependency of the solution. 

On the basis of the obtained results, we also investigated the secular increase in the 
astronomical unit (AU), reported by II17I : we showed that the leading order terms of the time- 
dependent component due to cosmological expansion are 9 orders of magnitude smaller than 
the observed value of dAU/dt, i.e., 15 ± 4 [m/cy]. Therefore, we explicitly asserted that it 
is not possible to explain the secular increase of AU in terms of cosmological expansion. 

Currently, it is assumed that the most plausible cause of dAU/dt is either the lack of 
calibrations of internal delays in the radio signals within spacecraft or the complexity in 
modeling the solar corona; however, no conclusive evidence has not been reported thus far. 
The origin of the secular increase in AU has also been examined in terms of other physical 
aspects such as solar mass loss I17II35I . the time variation of the gravitational constant G 
1171 , the influence of dark matter in the solar system [34 1, the multi-dimensional brane 
world scenario B6I . the transfer of rotational angular momentum of the Sun due to solar 
mass loss I37II . the azimuthally symmetric theory of gravitation l38l , and the kinematics in 
Randers-Finsler geometry 11391 . However, so far the cause of the secular increase in AU is 
not clear. Nonetheless, this phenomenon and other anomalies discovered in the solar system 
may be related to some fundamental properties of gravitation; therefore it is important to 
verify these phenomena theoretically and observationally. 

On the other hand, the cosmological effect on gravitationally bound local systems has 
been widely investigated, see i.e. 1301 and references therein. Its influence in the solar sys- 
tem is probably so small that one cannot detect it presently. Nevertheless, it seems that the 
theoretical discussion is not resolved. Therefore, from the theoretical point of view, it is in- 
teresting to construct robust theoretical model beyond the McVittie model and bring an end 
to the argument consistently. 

Owing to the rapid enhancement of observational techniques, more accurate and rigor- 
ous formulae are required to express light propagation. [10] and |11] succeeded in devel- 
oping an elegant and useful method for constructing the post-Newtonian/post-Minkowskian 
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approximation, and they showed that Synge's world function Q {xa,xb) (which is not widely 
adopted or discussed) is a powerful tool for formulations in practical problems. Thus, the 
value of Synge's world function can be reaffirmed, and it is possible to proceed with the 
theoretical development of observational models based on it. 
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